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ABSTRACT 

In a sample of elliptical galaxies that span a large range of mass, a previously unused Ca 
index, CaHK, shows that [Ca/Fe] and [Ca/Mg] systematically decrease with increasing elliptical 
galaxy mass. Metallicity mixtures, age effects, stellar chromospheric emission effects, and low- 
mass initial mass function (IMF) boost effects are ruled out as causes. A [Ca/Fe] range of less 
than 0.3 dex is sufficient to blanket all observations. Feature gradients within galaxies imply a 
global Ca deficit rather than a radius-dependent phenomenon. Some, but not all. Type II super- 
nova nucleosynthetic yield calculations indicate a decreasing Ca/Fe yield ratio in more massive 
supernovae, lending possible support to the hypothesis that more massive elliptical galaxies have 
an IMF that favors more massive stars. No Type II supernova nucleosynthetic yield calculations 
show significant leverage in the Ca/Fe ratio as a function of progenitor metallicity. Therefore, 
it seems unlikely that the Ca behavior can be explained as a built-in metallicity effect, and this 
argues against explanations that vary only the Type II to Type la supernova enrichment ratio. 



Subject headings: Galaxies: abundances - 
— Galaxies: stellar content 

INTRODUCTION 



Galaxies: elliptical and lenticular, cD — Galaxies: formation 



The determination of age and metallicity 
spreads of the stars that compose galaxies is highly 
useful for better understanding the formation and 
evolution of the content of the universe. Galax- 
ies that show little evidence for present-day star 
formation, such as elliptical galaxies, are better 
laboratories for unraveling the stellar populations 
than star-forming galaxies plagued by dust, neb- 
ulae, and bursty recent star formation. 

It is now well established that other galax- 
ies, large elliptical galaxies in particular, display 
chemical abundances differing that of our own 
Galaxy. Nonsolar abundance ratios in elliptical 
galaxies w e re firs t me ntioned as a pos sibility by 
lO'Connelll (Il976l) and iPetersonI (Il976l) . both of 
whom found Na and Mg features to be strong com- 
pared to Fe and Ca features. That [Mg/Fe] was 



truly greater than zero in large elliptical galaxies 
and not just a side effect of an overall lockstep en- 
hancement, w a s solidified by other s tudies (e.g., 
Peletieil Il989l: IWorthey et al.l Il992l: iFaber et al 
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2010I ). establishing the implication that there ex- 
ists an overall enhancement of some lighter ele- 
ments relative to iron. 

Magnesium and other even-Z elements lighter 
than the iron peak are made via capture of al- 
pha particles, starting with ^^C, and going on to 
include O, Ne, Mg, Si, S, Ar, Ca, and possibly 
Ti, where Ti might be too close to the Fe peak 
elements to be included. Workers tend to expect 
these elements to track one anoth er due to their 
comm on nucleosynthetic origin ([Wheeler et al. 
19891) . So Ca abundance should track Mg abun- 
dance according to this logic. 

Emerging observations of Ca indices, however, 
indicate that Ca tends not to follow Mg very 
well. Up to now, there have been four calcium 
feature indices used to measure calcium abun- 
dance in galaxies. These are Ca4227 and Ca4455 
(|Trager et all Il998l: IWorthevI Il998l ). the Ca II 
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"triplet" at » 8500 A (ICenarro et al.l l2003l). and 
Ca4227r (IProchaska et al.ll2005[ ). 



The Ca4227, Ca4455, and Ca II triplet indices 
suggest a mild decrease or a nea rly flat level of 



[Ca/Fe] with increasing [Mg/Fe] jVazdekis et al 
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19991: iSaglia et all 120021 : I Thomas et all l2003t 
Cenarro et al.ll2003l [20ollSmith et al.ll2009l ). In- 



terestingly and uniquely, the Ca4227r index pro- 
duces results to the contrary; higher obtained 
[Ca/Fe] abundances follow more clo sely the trend 
of [Mg/Fe] (Procha ska et al.ll2005[ ). That index 
was designed to cut out some CN absorption at 
the expense of requiring accurate relative spec- 
trophotometry. 

In terms of observed indices versus galaxy 
velocity dispersion, Mg-sensitive indices exhibit 
strong correlations with galaxy velocity disper- 
sion (thus, roughly, galaxy mass), while a much 
more modest co rrelation exists for Fe-peak indices 
()Worthevlll998h . The Ca II triplet shows a mod- 
est anticorr elation between C a and central velocity 
dispe rsion (jSaglia et al.ll2002l : ICenarro et al.ll2003 . 
20041 ) ■ whil e the trend with Ca4227 and Ca4455 is 
nearly flat (iTrager et al.lll998l:[Cenarro et al.l2004 : 



Smith et al JI2OO9I) Sample size and velocity disper 
sion range from 28 to 147 galaxies and a ~ 40—370 
km s""'^, respectively, with no obvious trend in 
results de pendent on either of these sample pa- 
rameters (ISas . lia et al. 20o3 Thomas et al. 2003 : 



Cenarro et al.ll2003i 120041: ISmith et al.ll2009( ). The 
trend of Ca4227r i s increasing with inc reasing 
velocity dispersion (jProchaska et al.l 120051 ) . how- 
ever, the authors limited their sample of 175 
early-type galaxies to those with a velocity dis- 
persion of cr < 230 km s~^ out of the full range of 
cr = 50 - 300 km s^^. 



An explanation discussed bylSaglia et al.l (|2002l ) 
and lCenarro et al.l (|2003ll2004l ) involves the steep- 
ening of the initial mass function (IMF) at the 
low-mass end to ma ke more dwarfs. This w ould 
increase Na strength ( Spinrad fc Tavloilll97l[ ) and 
decrease CaT strength, both of which phenom- 
ena are observed in high-mass elliptical galax- 
ies. While it seems appealing to solve two ob- 
servational effects with one astrophysical reason, 
boosting the number of low-mass stars also greatly 
increases the near-IR flux, reddening colors involv- 
ing near-IR passbands, and lowers near-IR M/L 
ratios, and drastic lower-IMF changes are usu- 



ally looked upon with ill favor on tho se grounds 
(jSaglia et al.ll2002HThomas et al.ll2003[ ). 

A more mainstream explanation rests on the as- 
trophysics of supernovae. Generally, Type la su- 
pernovae produce mainly Fe-peak elements, while 
Type II supernovae produce the lighter metals, 
and alpha capture elements in p articular. What- 
ever the causes [many ar e listed in lWorthevI ( 1998h 
and iTrager et al.l (|l998t ) and include time delays, 
wind scenarios, IMF changes at the high-mass end, 
and binarism trends], the Fe-peak could plausibly 
decouple from the light elements in terms of abun- 
dance, but the light elements should be in lockstep 
with each other. Thus, the separation of the be- 
havior of Mg and Ca is inherently puzzling. 

Worthev (Il998h proposed that the simplest 



explanation for the Ca flatness might be due 
to metallicity-dependent supernova yields, where 
the production or ejection of Ca is decreased in 
metal-rich superno vae. This is weakly suppo rted 
by the studies of Wooslev fc Weaven ( 1995 ). in 
that [a/Ca] increases slightly with increasing 
metallicity, but more recent calculations from 



Nomoto et al.l (|2006f ) show no evidence of this 
trend, with Ca tracking Mg in supernova yields of 
all metallicities. The latter study includes metal- 
licities greater than solar. With Occam's razor 
thus tentatively defeated, some other effect must 
be operating to create a trend along the mass se- 
quence of elliptical galaxies, such that low-mass 
ellipticals, lenticulars, and spirals have a nearly- 
solar mixture, while very large ellipticals (that 
sit alone at the high-mass end of all galaxy mor- 
phological types) have a light-element enhanced 
mixture. 

Both the observational material and the mod- 
els used vary quite widely among the workers that 
have addressed the Ca topic. Since the conclu- 
sions also vary, this motivates the measurement of 
[Ca/Fe] in a better way, if possible. In this paper, 
we approach the calcium puzzle with the Ca II 
Fraunhofer H and K features that are strong fea- 
tures in cool stars. Elliptical galaxies are largely 
composed of low-mass, cooler stars, and thus the 
integrated-light spectra of these galaxies also dis- 
play these strong features, leading to a nearly ideal 
index for measuring calcium abundance trends 
if the stellar features are indeed sensitive to Ca 
abundance. Although Ca H and K are subject 
to absorption by the interstellar medium, ellipti- 
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cal galaxies contain very little interstellar matter. 
The blue wavelength of the feature ensures that 
ordinary G and K type stars, mostly main se- 
quence and warmer giants, contribute most of the 
light, in contrast to the CaT 8600 feature, which 
is dominated by red giant branch (RGB) tip and 
asymoptotic giant branch (AGB) star light whose 
numbers, temperatures, and s pectra are all l ess 
well understood (JBuzzoni 1989i: IWorthev 1994 ). 

In this paper we focus on a partic ular CaH+K 
index, defined in IServen et al. I (I2OO5I) and named 
CaHK, and discover for the first time an emphatic 
anticorrelation between CaHK and Mg b index 
strength. We follow in §3 with an exploration of 
the influence of calcium emission lines on CaHK 
index strengths and rule out Ca H emission rever- 
sals as a possible culprit for calcium underabun- 
dance trends. Admixtures of metal poor subpopu- 
lations or young subpopulations or dwarf-enriched 
populations are also ruled out, leaving abundance 
ratio changes as the only viable explanation. In 
§4, we add galaxy-internal gradient information 
and discuss the inferred abundance ratio trends 
between galaxies and within galaxies. Finally, §5 
summarizes the findings and discusses them within 
the context of basic chemical evolution. 

2. MODELS AND DATA 

A grid of synthetic stellar spectra resampled 
to O.5A wavelength intervals in the ranM 
to lOOOOA was incorporated into the IWorthevI 
( I994I ) infrastructure (in which evolution choices 
and spectra l library choices are modular ; cf 



Dotter et all 120071 : iLee et al.l 120091 : IServen et al" 



20101 ) in order to produce integrated-light mod- 
els. This spectral grid was repeatedly calculated 
with small abundance changes, for example, a 
0.3-dex increase in calcium abundance, relative to 
the scaled-solar mixture, so that spectral shape 
changes due to elemental abundance drifts could 
be calib rated. A list of 23 e l emen ts was included, 
and the Grevesse &: Sauval (Il998l) abundance list 
was used as the solar mixture. 

A pertinent example of this is given in Fig. 
[1] which shows the Ca-enhanced spectrum di- 
vided by the unenhanced spectrum. Several Ca- 
sensitive features become obvious. Ca does not 
form molecules in stars, so the features are due 
to Ca I or Ca II atomic species. It is obvious 



from Fig. [T] that the Ca H and K features give 
the strongest percentage response to increased Ca 
abundance. Furthermore, the response is spread 
over a quite broad wavelength regime compared 
to the other Ca-sensitive regions of t he spectrum 



making measurements even easier. IServen et al 
POOSI ) noted that CaHK is a factor of three more 
sensitive than the CaT indices despite (in their 
formalism) being at much lower signal to noise ra- 
tio (S/N). For the same S/N, CaHK is almost an 
order of magnitude better. 

In order to make an integrated light model pop- 
ulation spectrum, new stellar index fitting func- 



tions were generated (first cited in iPoole et al 



2010h . The three sources of stellar spectr a l data 
were (1) a variant of the IWorthev et al. (1994) 
Lick spectra with a wavelength scale refined via 
cr oss-correlation, (2) the MILES spe ctral library 
of Sanchez-Blazauez. P. et all (l2006bl) . and (3) the 
Indo-US librarv of lValdes et al.l t004) . The Indo- 
US library, smoothed to mimic a Gaussian velocity 
dispersion with a — 200 km s~^, was adopted as 
the fiducial set, with the other two adjusted to 
fit via linear transformations. Both sets, in gen- 
eral, needed more than merely zero-point shifts. 
In the Lick slope term was sometimes used 

to account for the difference in spectra resolution. 
In the MILES case, a slope term was often ap- 
plied to account for what is probably a scattered 
light effect, which is discussed in g reater depth in 



Sanchez-Blazauez. P. et al.l ()2006 



The data were then fit with multivariate poly- 
nomials in five generously overlapping tempera- 
ture swaths, then summarized in a lookup table 
that gives each index strength as a function of 
stellar atmosphere parameters. As described so 
far, the fits mimic only the run of Galactic star 
index strengths with a single abundance parame- 
ter, [Fe/H]. In order to add the sensitivity to indi- 
vidual elements, the synthetic spectra were again 
employed to generate index responses dl /d[X/K], 
where X is the element of interest, mostly calcium 
in this paper, and R stands for "generic heavy el- 
ement," and is employed so that the common case 
of Fe abundance changes can be written as [Fe/R] 
and will still make logical sense. 

This library of element-sensitive indices was 
then attache d to stellar evo lutionary isochrones 
by choosing a lSalpeten (|l955l ) IMF. The isochrone 
sets used in the modules are modular, but initial 
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Fig. 1. — Ratio of two model spectra of simple 
stellar populations of age 15 Gyr and abundance 
around solar, smoothed somewhat for clarity. The 
model with calcium raised a factor of two is di- 
vided by the unenhanced model. The area of the 
spectrum that responds best to extra Ca is the Ca 
HK feature, decreasing 30% over a broad wave- 
length range. Lesser spectral responses at the red 
calcium "triplet" and at 4227A are noted. 



testing was done with the lWorthevI (|l994l ) set. To 
produce an integrated-light index strength, values 
of the flux and continuum for each "star" in the 
isochrone, weighted by number, are added up the 
isochrone, and then reformed into an index value 
after summation. Finally, since we compare in the 
end to spectra that are smoothed to mimic a 300 
km s~^, index corrections were calculated using 
synthetic model spectra smoothed to the relevant 
velocity dispersions and applied to the model in- 
dices. 

For observational comparison, we use four data 
sources for the inde x measurements. They are the 
Graves et ahl (|2007t ) SDSS averages of non-LINER 
elliptical galaxies and three sets of high-S/N near- 
nuclear spect r a (ISanchez-Blazguez et al.l 2006at 



Serven et al.l l20ld iTrager et al 



Serven 120101 : 
2005, ,20081 ). The Graves average spectra are 
parceled out into bins of velocity dispersion, cov- 
ering a wide r ange. The CaHK index is defined in 



are defined in 



Serven et al.l ( 20051) and the regu lar Lick indices 



Trager et al.l (|l998[) . 



The iTrager et al. ( 20081 ) spectra were kindly 



provided by the lead author of that paper, already 
smoothed to 300 km s~^. They are 12 early- type 
galaxies in the Coma cluster with velocity disper- 
sions between 41 and 270 km s~^, including the cD 
galaxy NGC 4874, taken with the Low-Resolution 
Imaging Spectrometer on the Keck II 10 m Tele- 
scope. This instrument uses slitlets as spectral 
apertures, and the reduced spectra were weighted 
along the slit to mimic a circular aperture of diam- 
eter 2.7". The S/N in a 75 km s""'^ bin was between 
35 and 292, and one fiux standard was observed for 
purposes of applying a spectral shape correction. 
W e were also kindly provid e d a sub set of 35 of the 
98 ISanchez-Blazguez et al. I (l2006ah long slit spec- 
tra, trimmed to a 4"equivalent (at redshift z = 
0.016 for a physical size of 0.62 kpc assuming Hq — 
70 km s^^ Mpc^^) central aperture. They were 
obtained at the 4.2m William Herschel Telescope 
at Roque de los Muchachos Observatory using the 
ISIS spectrograph and the Cassegrain Twin Spec- 
trograph on the 3.5m German-Spanish Astrophys- 
ical Observatory. The 35 elliptical galaxies cover a 
range of velocity dispersion between 130 and 330 
km s~^ and also local density. The spectra are 
at a range of native instrumental resolution be- 
tween 3.5A and 6.6 A (FWHM) and a typical S/N 
of 50-110 A^^. A few flux standards were observed 
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during the runs for purposes of spectra shape cor- 
rection. We note that CaHK is near the bhic edge 
of the spectral cove rage. The Serven ^rvcn 2010!; 
Serven et ah 201Cl[ ) spectra of mostly Virgo ellip- 
tical galaxies span a range of velocity dispersion 
from 80 to 360 km s^^ and were obtained with the 
Cassegrain Spectrograph at the 4 m Mayall tele- 
scope at Kitt Peak National Observatory. Nuclear 
portions of the long slit spectra were extracted at 
an aperture of 13.8", corresponding to a physical 
distance of 1.1 kpc assuming a Virgo clust er dis- 
tance of 16.47 Mpc (Blakcsl ce et al.ll2009t) . The 
S/N is over 100 per 2A bin except for the dwarf 
galaxies in the sample. Dozens of flux standard 
star observations of several flux standards at vary- 
ing airmasses and with both narrow and wide slits 
were observed nightly during the KPNO runs for 
spe ctra shape cor rection. Of the various data sets, 
the Serven ( 2010f ) spectra are most likely to have 
the most reliable spectrophotometric shape, due to 
the many flux standards observed and also because 
the CaHK index falls in the middle of the spectral 
range of the instrument where scattered light was 
demonstrably negligible, even for red objects, and 
camera focus issues were never observed. 

The apertures sample beyond the half-light ra- 
dius of each SDSS galaxy but are more nearly nu- 
clear for the high-S/N sets of local galaxies. For 
galaxies (or average galaxies) with velocity disper- 
sion greater than 300 km s~^, we computed correc- 
tions via model spectra the same as for the model 
indices described above but in the opposite direc- 
tion. For galaxies with velocity dispersion less 
than 300 km s~^, extra, wavelength-dependent 
Gaussian smoothing was applied to bring them to 
300 km s~^ before indices were measured. 

The CaHK index definition is quite broad to 
capture most of the absorption signal that is 
present. This broad wavelength span leads to 
probable zero-point offsets between data sets due 
to imperfections in recovering relati ve spectropho- 
tomet ric shape, as discussed in .Worthev fc Ottaviani 
(|l997l ). There a r e only two galaxies in common 



ily 



between Serven (20101) and the subset of the 



Sanchez-Blazquez et al.l (12006a|) with which we 
worked, and within the possible offsets between 
the sets due to aperture effects and small num- 
ber statistics no mutual zero-point offset in the 
CaHK index can be unambiguously shown. No 
other data sets have overlap. 



Some basic index behaviors against velocity dis- 
persion are shown in Figure O The Mg b index 
correlates strongly with velocity dispersion, with 
wide range (2A). Indices <Fe>= 0.5 x (Fe5270 
-t- Fe5335) and Ca4227 correlate at lesser signifi- 
cance. The Ca HK index, in contrast, anticorre- 
lates. The offsets between data sets, highlighted in 
the bottom panel of Figure [2] by splitting the data 
into the four separate sets and fitting lines sep- 
arately, are predominantly due to (relative) spec- 
trophotometric shape effect s. The data are fit with 
the F/TE'Xr program from lPress et all ()l992l ) and 
the resulting two-error least-squares fits are sum- 
marized in Table [T] 

The behaviors are as observed in the past for 
the three top indices, all positive correlations, with 
Mg b being the strongest correlation and Ca4227 
the weakest. The CaHK slope, however, is neg- 
ative in all four data sets. The zero point shifts 
we take in stride for this work since our conclu- 
sions are not affected by them, or even by mod- 
est slope differences, and we have not tried to 
aperture-correct the sets. (Such corrections would 
be stronger for Mg b and <Fe> and weaker for 
CaHK, based on observed gradient strengths.) 

With these models and data in hand, we dis- 
play them together in Figs. [3]and|4l arraying the 
CaHK index against the <Fe> index in Fig. [3] 
and against the Mg b index in Fig. SI A second 
panel is provided in each case, with deltas applied 
to the data to make the four data sets agree with 
the average value at a velocity dispersion of 250 
k m s~^. As in quite a few such index-index plots 



(jWorthevI 119981 ) . the galaxies follow a trajectory 
significantly skewed from Galactic expectations, 
indicating systematic drifts in abundance ratios 
away from Galactic trends. In this case, the sense 
is that the stronger lined elliptical galaxies have a 
weaker calcium index. Holding age fixed, stronger 
metallic lines imply higher metallicity. Further- 
more, in the case of Mg b, Mg b correlates wit h 
velocity dispersion strongly ([Bender et al.lll993l) . 
implying that it is the larger galaxies that have 
the weakest CaHK index. This can be seen via 
the symbol size coding in Figs. [3] and |4] well; ex- 
cept, perhaps, for the Herschel data set, the small 
galaxies lie near the models while the large ones 
are far away. 

Iron indices have a weaker correlation with ve- 
locity dispersion, indicating that larger ellipticals 
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Table 1 

Fits Versus Velocity Dispersion 



Diagram 


Intercept (A) 


Slope (A per km s ^) 


Mg 5, combined 


2.55 ±0.07 


0.0070 ± 0.0003 


<Fe>, combined 


1.79 ±0.04 


0.0023 ± 0.0002 


Ca4227, combined 


0.92 ±0.04 


0.0010 ±0.0002 


CaHK, KPNO 


19.99 ±0.32 


-0.0226 ±0.0018 


CaHK, SDSS 


22.19 ±0.43 


-0.0091 ± 0.0025 


CaHK, Keck 


20.83 ±0.24 


-0.0123 ±0.0012 


CaHK, Herschel 


25.90 ±0.12 


-0.0241 ± 0.0006 




100 200 300 

Velocity Dispersion (km s"') 



Fig. 2. — Four feature strength indices vs. 
galaxy centra l velocity dis p ersion . Data from 
Serven (2010*); Scrvcn et al. (20ld) (red squares), 
Sancl icz-Blazqucz et al. (2006a) (light blue open 
circles), ICraves et all (|2007t ) (green stars), and 
Trager et all (| 2008>) (filled blue circles) are plot- 
ted with symbol size proportional to velocity dis- 
persion. The lines are two error linear fits to the 
combined data, except in the bottom panel, where 
four lines are shown, split by data set to show sys- 
tematic offsets. 



have only slightly more Fe than small ones. That 
trend is weak, and, indeed, it is still in the realm of 
possibility that age systematics might be able to 
drive the Fe trend. Youth drives the stellar popu- 
lation left and down in both of these plots, as seen 
in the model grid lines. 

3. EXPLANATIONS FOR CALCIUM 
INDEX BEHAVIOR 

From the models, it is clear that gross age 
and metallicity effects cannot account for the div- 
ing CaHK index since mixtures of various simple 
stellar population add approximately like vectors 
within the model grid, but the data points lie out- 
side the grid, even if one were to shift the grid 
artificially. This leaves four possibilities. One, 
chromospheric emission (jWilson Sz Vainu Bappu 
Il957h may operate with sufficient strength to be 
an issue. Two, the wavelength difference between 
CaHK around 3900A and the 5100A to 5400A re- 
gion where Mg b and the Fe indices lie may leave 
the door open for a small fraction of metal-poor 
subpopulation to dilute the CaHK index, with the 
more metal-rich, larger galaxies having less blue 
fiux and hence more fill-in with the weak-lined, 
blue metal-poor subpopulation. Three, an en- 
riched dwarf population within the galaxy creates 
a shallower IMF on the low-mass end and weakens 
the CaHK index as it might for the Ca II triplet 
feature. Four, the Ca abundance relative to Fe 
and Mg is changing, diluting with increased galaxy 
size, with implications for supernova enrichment of 
the galaxies. 
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3.1. Chromospheric Emission 





Fig. 3. — CaHK is sliown a s a func ti on o f 
<Fe>. Models based on tlie IWorthevI (|l994l ) 
isoclirones are sliown, at ages 8, 12, and 17 
Gyr, witli isometallic lines drawn at scaled 
solar [M/H] every 0.25 dex. Observational 
points are drawn with symbol size propor- 
tional to_galaxy_veloci^^ Data from 
ServenI (|2010|) : IServen (120101) (red squares), 
Sanchez-Blazguez et al. I (l2006al) (light blue open 
circles). Graves et al.l ()2007[ ) (green stars), and 
Trager et all (|2008[ ) (filled blue circles) are plot- 
ted. In panel (a) the data are as measured. In 
panel (b) the CaHK values have been shifted to 
the average value at 250 km s^^ from the index 
vs. velocity dispersion relations of Table 1. 

















1 — ^ 





















\ 




Mg b (A) 



Fig. 4. — CaHK is shown as a function of Mg h. 
Symbols and lines are as in Fig. [3l 



Within the centers of the H and K lines, 
emission spikes are often present for stars of 
ty pe GO and later. This was fir st discovered 

by ISchwarzschild fc Eberhard ( 19131) and was fol- 

lowed up by extensiv e literature. Wilson fc Vainu Bappul 
(|l957l ) and IWilsonI (ji959) found the logarithms 
of the widths of these emission lines to be lin- 
early correlated with luminosity, independent of 
spectral type, although stars of a given luminos- 
ity could have a wide range of intensities of the 
emission. This emission desc ribes general chromo- 
spher ic activity in stars (iWilson &: Vainu Bappu 
19571 ) and is of interest because it results in weaker 
Ca H and K equivalent width index measurements. 
We cannot use the previous work directly because 
we need to judge its effects on the integrated light 
of stellar populations in the particular index that 
we are using. 

Therefore, we proceeded to estimate the max- 
imum amount of emission and its effect on the 
CaHK index by examining a n observationa l spec- 



(|2004l ) 



tra library, choosing that of iValdes et al. 
due to spectral resolution and wavelength cover- 
age. For purposes of this exercise, the effective 
temperatures of the stars were adopted right from 
the library. Giants of various temperatures, and 
excluding metal-poor or chemically peculiar stars, 
were visually examined for emission in the cores of 
the lines. If emission was spotted, a hand-drawn 
guess for the true depth of the core was inserted 
into the spectrum, and the CaHK index measured 
both with and without the emission spike. 

Figure [5] shows the spectra, before and after, 
of a star with strong emission, HD 200527. Af- 
ter the star is "corrected," it is smoothed to a 
galaxy-like velocity dispersion (200 km s~^) that 
serves as the standard for the index system. Fig- 
ure [6] shows the results for a collection of K and 
M giants, specifically ACaHK — CaHK (corrected) 
— CaHK (uncorrected). While, as expected, there 
is quite a lot of scatter in the intensity of the 
emission, there is also a fairly well-defined max- 
imum edge. We examined many stars hotter than 
4500 K, but none of them had any convincing Call 
emission. The late-type giants examined here have 
CaHK index values of order 30A . 

In order to complete the test, we approximated 
the ACaHK correction to be a line that is 0.5 A at 
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Fig. 5. — Representative set of spectra showing Ca 
H (A 3969) and K (A 3934) features. The plot hm- 
its are precisely the index passband of the CaHK 
index definition, while the pseudocontinua pass- 
bands are outside the limits. The solid line is the 



spectrum from the Indo-US library of lValdes et al 



|2004) for M giant HD 200527. The dotted por- 
tions of the spectrum show the altered Ca H and 
K lines to remove emission. The dashed line is the 
latter spectrum smoothed to mimic a Gaussian ve- 
locity dispersion with a — 200 km s^^. 



3000 K and 0.0 A at 4500 K and hotter, and then 
a pplied this functi on within the population models 
of lWorthevI (|l994l ). At ages 3, 6, and 12 Gyr, the 
AGaHK values were 0.015A , 0.016A , and 0.024A 
, respectively. The small magnitude of these shifts 
results from the cool onset of the Call emission 
phenomenon, and the rapidly falling continuum 
levels in these stars, so that these stars have very 
minimal contribution at the blue wavelengths of 
Ca H and K. Definitively, then, the hypothesis 
that Call emission has anything to do with the 
peculiar morphology of Figs. [3] and |4] is ruled out. 

3.2. Metal-Poor Subpopulations 

The second hypothesis, that a varying ratio (in 
light at 4000 A) of metal-poor to metal-rich sub- 
populations can account for the galaxy trend can 
also be rejected. Mixtures of populations add ap- 
proximately like vectors within the model grid. 
Therefore, one is generally trapped within the con- 
fines of the model grid when mixing different pop- 
ulations together, the grids in Figs. [3] and |4] ex- 
tend between [Fe/H] = -2 and [Fe/H] = -1-0.5 
and lie on a roughly linear, roughly orthogonal 
trajectory. There is no way to produce the galaxy 
trend with inventive metallicity mixtures, with the 
caveat that the fairly large wavelength separation 
between Mg h and <Fe> in the green and CaHK in 
the violet might possibly introduce a nonlinearity. 

To investigate this, we concoct some more com- 
plex models in Fig. [T] The sequence marked 0%, 
5%, and 10% in Fig. [7] is a two-metallicity com- 
posite stellar population with the mass fraction 
of the metal-poor component ([M/H] = —1.5) as 
labeled. The age is 8 Gyr for both components, 
and the metal-rich component starts at [M/H] = 
and extends to [M/H] = 0.25 at the end of the 
line segment marked with the letter "Z." The phe- 
nomenon of note is the slight twisting of the three 
line seg ments as the metal-p oor fraction increases. 



See also Serven et al 



(l2010l) . who explore indices 
that span even further in wavelength. 

The twisting is clockwise, and that is in the 
correct sense to make the model sequences fiat- 
ten or even decrease with increasing metallicity. 
The physical reason is that metal-poor popula- 
tions, dominated by the main sequence stars at 
this wavelength, are weak-lined: as the metal- 
rich component drifts yet more metal-rich, its 
fractional light contribution in the blue decreases 
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dramatically, giving more weight to the metal- 
poor component. There are two reasons, how- 
ever, not to take this very seriously. First, the 
model would imply that more metal-rich galaxies 
have a more broad abundance distribution. There 
is no evidence for th is at less than solar [M/H] 
(|Worthev et al.ll20n,'^^ . and, one would think, little 
motivation to demand such a behavior at higher 
abundance. Second, metal-poor fractions as high 
as 10% are probably ruled out by other observa- 



tions with less than 5% preferred (jWorthev et al 



[1996) and this leaves the model twists at a very 
mild level. 

Figure [7] also shows the model vectors for bulk 
age increase and lockstep abundance increase, for 
basic orientation. It also shows three vectors for 
Ca, Fe, and Mg abundance mixture tweaks that 
are useful for the discussion that follows. 

3.3. Dwarf Enrichment 



Th e boo sted lower IMF hypot hesis of lSaglia et al, 
([2002) and lCenarro et al.l (j2003l ) can be addressed 
with the CaHK index as well. Models that include 
lower main-sequence stars down to 0.1 M©, of age 
12 Gyr and of solar composi tion were computed. 
These models used lWorthevI ( 19941 ) isochrones due 
to the availability of low-mass stars in them. The 
IMF was assumed to be a power law. Slopes of 
2.35 (ISalDeteijll95"5h and 5.35 produced CaHK in- 
dex values of 22.48A and 21.04A , respectively. 
However, for this ACaHK ~ 1.5A, the same two 
IMF slopes gave V — K colors of 3.45 mag and 4.99 
mag, respectively. Dust free elliptical galaxies do 
not have colors much redder than V — K ss 3.6 
mag. To get a ACaHK from -0.5A up an or- 
der of magnitude to the ^5A needed to match 
the observations would require IMF slopes much 
greater than 5.35. That seems impossible, given 
the strong infrared color constraint, even if the 
lower-mass end of the IMF was tapered and fine- 
tuned to decrease the number of cool dwarfs. 



Th e more modest proposal of lvan Dokkum fc Conrov 
(|2010 ). who consider slopes of 2.35, 3.0, and 3.5 
with the same 0.1 Mq lower cutoff and compare 
to Wing-Ford and Na I features in the red would 
probably not violate photometric constraints and 
would be consistent with our observations. We 
note that those authors did not consider the hy- 
pothesis of enhanced Na abundance, correction 
for which would drive their results toward more 



Salpeter-like slopes. 

3.4. Element Ratio Trend 

The explanation that remains is that the Ca 
abundance is changing systematically from smaller 
elliptical galaxies to larger ones . This is simila r to 
the trends for Mg, Na, and N, (|Worthevlll9"98l ) al- 
though the Ca abundance decreases instead of in- 
creases. These trends clearly arise through chem- 
ical enrichment processes, and primarily through 
supernova yields. 

The magnitude of the [Ca/R] abundance spread 
depends on how conservative one would like to 
be when interpreting Figs. |3l IH and [T] If the 
model-observation separation and >10A CaHK in- 
dex spreads are to be taken at face value, then the 
[Ca/R] variation must be at least —0.15 dex for 
the weakest-lined galaxies, as illustrated by the 
vertically-down Ca vector in Fig. [T] Naturally, 
the trends in [Fe/R] and [Mg/R] must also be 
taken into account, and those elements do affect 
the CaHK strength, by almost 3 A for a [Mg/R] 
shift of 0.2, sufficient to explain the Mg-cr trend. 
At face value, that is not sufficient to avoid in- 
voking a [Ca/R] change. (Note that a positive 
A[Fe/R] implies a negative A[Ca/Fe].) 

If one, not unreasonably, would rather put zero 
faith in the observational and model zero points, 
then one is left with relative change vectors, and 
only the massively-average SDSS data set should 
be compared with relative model change vectors. 
This is enabled in Fig. [71 where, looking only 
at the SDSS data and giving complete zero-point 
freedom to the models, the vertical change can be 
explained with the 0.2 Mg enhancement, with the 
need for an [Fe/R] boost of about twice what is 
illustrated; about 0.1 dex in [Fe/R] (or [Ca/Fe] 
= —0.1) while leaving [Ca/R] unchanged. Or, 
one could effect the same change with depleted 
Ca. Either way, and under liberal or conservative 
assumptions, [Ca/Fe] does decrease from small 
galaxies to large. 

As regards calcium, expectations from nucle- 
osynthesis theory are contradictor y at the mo- 
ment. If we look at the work of Nomoto et al 



(|2006l ) and iKobavashi et al.l (|2006l) we find that 
there is a strong metallicity dependence in Na 
yield that has not been predicted before. This 
could easily explain the high Na abundance in 
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massive elliptical galaxies. However, the Ca yields 
track Ti and Mg very well (normalized by Fe 
yield), both as a function of supernova progeni- 
tor mass and as a function of supernova progenitor 
initial abundance. However, if we look at the work 



of IWooslev &: Weaveij (|l995f ) and iTimmes et al 



(|l995l ) we find a very strong dependence of Ca 
yield with supernova progenitor mass: Ca de- 
creases relative to Fe, while Ti is about constant 
and Mg increases with increasing progenitor mass, 
and the variation is more than the required fac- 
tor of two. This is in the correct sense to ex- 
plain our observations if IMF variation (at the 
high mass end) is the root cause of the abundance 
trend. In this scenario, a galaxy with higher ve- 
locity dispersion will provide star formation en- 
vironments that are biased toward higher mass 
star formation. This will chemically enrich the 
g alaxy with the highe r mass progenitor yields, and 
if IWooslev fc Weaveij ([l99i) yields are used, then 
it could explain the trend we see. 

The other favored explanation for the observed 
[Mg/Fe] trend is a timescale argument in which 
a longer timescale for enrichment causes Type la 
supernovae to contribute substantially, adding Fe- 
peak elements to the mixture. However, Type 
la supernovae do not add enough Ca or Ti or 
lighter elements to generate observable abundance 
trends within the lighter metals. Therefore, the 
Ca trend is not readily explainable in this hypoth- 
esis, because no published, detailed set of super- 
nova yields has sufficient metallicity dependence 
(initial metallicity of the progenitor stars, that is) 
to torque the [Ca/Mg] ratios into what we would 
infer. The Ca evidence therefore tends to disfavor 
the timescale argument and favor the IMF argu- 
ment. 



cases. Therefore we predict that the Ca dilution is 
less than a factor of two relative to a scaled-solar 
abundance mixture. 

Figs. H and El also show gradient trends from 
the IServenI (|2010|) spectra, where the near- nuclear 
regions are averaged and plotted as crosses, and 
the more distant regions are averaged and plot- 
ted as the endpoints of trend lines. The endpoints 
are the most distant points available with good 
data, and the end distance is quite variable, but 
the main point of the figure is to display the most 
reliable vector direction. With some exceptions, 
the gradients tend to lie along the age-metallicity 
direction indicated by the model grid, and not 
along a direction more indicative of strong internal 
gradients in the [Ca/Fe] or [Ca/Mg] abundances. 
This is in agre ement with previous observations of 
Mg gradients (jWorthev et a l. 19921 and indicates 
a rather global enrichme nt pattern through o ut the 



gala xy. More recently , Kuntschner et al. ( 2006[ ) 



and iRawle et al.l (|201Cl| ) present observational and 
model-dependent results, respectively, that de- 
scribe steady, or, on average, perhaps slightly neg- 
ative Ca abundance gradient trends with increas- 
ing galaxy radius. The data presented here are 
certainly not sufiicient to rule out those conclu- 
sions. 

A straightforward way to homogeneously enrich 
the volume within a galaxy is to have fairly effec- 
tive spatial mixing during the enrichment process. 
Since the process itself is almost certainly super- 
nova blast waves, this seems reasonable. Other 
mechanisms could also contribute, such as phase 
mixing of the stars by dynamical effects long after 
the enrichment phase is over. 

5. SUMMARY AND CONCLUSION 



4. INTERIOR GRADIENTS 

One of our data sets contains spatially resolved 
indices. In Figs. [5] and [3] we show an additional 
model grid that has been altered such that Mg is 
enhanced by a 0.3 dex, and Ca is decreased by 0.3 
dex, so that [Ca/Fc] = -0.3 and [Ca/Mg] = -0.6 
and [Mg/Fe] — +0.3. The heavy element mass 
fraction is not (quite) preserved in this altered 
mixture, and increases a bit because Mg is more 
abundant than Ca. This modest amount of alter- 
ation is sufiicient to explain even the most extreme 



Four samples of elliptical galaxies that span a 
large range of mass are compared to new models 
that feature a previously unused Ca index, CaHK, 
and that are also quite flexible in displaying chem- 
ical abundance effects on spectra and spectral in- 
dices for integrated light. It is seen that [Ca/Fe] 
and [Ca/Mg] systematically decrease with increas- 
ing elliptical galaxy Mg and Fe line strengths, and, 
by implication, galaxy mass. Metallicity mixtures 
and age efi'ects are ruled out as causes due to in- 
ability to match the observations. Stellar chro- 
mospheric emission was explored but rejected as 
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being about two orders of magnitude too weak to 
have any effect on the observed trend. Modula- 
tions of the lower IMF increasing the dwarf light 
fraction do not affect the CaHK index, and thus 
low-mass IMF variations are also ruled out as a 
cause. 

That leaves chemical abundance variation as 
the sole remaining explanation. The size of the 
abundance shifts are modest. Less than a factor 
of two in Ca dilution can easily explain the trend. 

Feature gradients within galaxies imply a ten- 
dency toward a pangalactic Ca deficit rather than 
a strongly radius-dependent phenomenon. The el- 
ement mixture trend may match a predicted su- 
pernova progenitor mass dependence on the Ca 
yield, although this prediction is not yet firm, ap- 
parently. If the mass-dependent Ca yield is con- 
firmed, then the observations strongly support a 
variable IMF (at the upper end only) where more 
massive galaxies tend to form more of the most 
massive stars. In any case, the Ca trend does not 
particularly favor the scenario of timescale modu- 
lated Type II vs. Type la supernova enrichment 
effectiveness that has received most mention in the 
literature to date. 
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Science Foundation grant AST-0346347. The au- 
thors thank E. Toloba, P. Sanchcz-Blazqucz, and 
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Fig. 6.— Trend of ACaHK = CaHK (corrected) 
— CaHK (uncorrected) vs. Toff for cool metal- rich 
giant stars selected from the Indo-US library of 
Valdes et all (l2004h . 



This 2-column preprint was prepared with the AAS lAT^rjX 
macros v5.2. 
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Fig. 7. — Modeling results in the CaHK vs. 
<Fe> (panel (a)) and Mg b (panel (b)) diagrams 
are shown. The shifted versions of the data are 
repeated from previous figures, sans error bars. 
Three heavy dots mark models of single-age pop- 
ulations, that age being 8 Gyr. The 0%, 5%, 
and 10% sequence refers to the fraction by mass 
of an additional metal-poor component of [M/H] 
= —1.5. The lighter dots indicate tweaks from 
that initial starting model. The ones marked "Z" 
refer to increasing lockstep heavy element abun- 
dance by 0.25 dex. In the case of the metallicity- 
composite models, only the metal-rich, dominant 
component changes its abundance. The other la- 
bels have the following meanings. "Age" means 
an age increase from 8 Gyr to 12 Gyr, "Frosting" 
is a 10% by mass addition of a 1 Gyr age, solar 
abundance component, "Ca— " is a drop of [Ca/R] 
of -0.15, "Fe-I-" is an increase of [Fe/R] of +0.05, 
and "Mg-K" is an increase of [Mg/R] of +0.2. 
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Fig. 8. — CaHK and < Fe> gradient trends. Using 
only the IServenI t01(t data set, gradient trends 
within each galaxy are plotted with line segments 
connecting averages of the near-nuclear spectra 
with averages of the outer parts of the light along 
the slit. The model grids contain ages 1.5, 2, 3, 5, 
8, 12, 17 Gyr populations and metallicities [M/H] 
= -2, -1.5, -1.0, -0.5, -0.25, 0, 0.25, and 0.5, 
with no ages less than 8 Gyr plotted for metal- 
licites less than —0.25. The lower set of models is 
identical except that, at the spectral level, 0.3 dex 
more Mg and 0.3 less Ca is included. This alters Z 
in a slightly positive way, but leaves all heavy el- 
ement ratios unaltered except those involving Mg 
or Ca. 
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Fig. 9. — CaHK and M r b gradient trends. Using 
only the IServenI t()l(i \ data set, gradient trends 
within each galaxy are plotted with line segments 
connecting averages of the near-nuclear spectra 
with averages of the last reliable outer measure- 
ment along the slit. Models are as in Fig. [H NGC 
3156 is dominated by a ~1 Gyr young stellar popu- 
lation and also appears to buck the gradient trend 
defined by the other galaxies somewhat. 
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